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Abstract 



We theoretically studied the electronic and electrical properties of metallic 
and semiconducting peapods with encapsulated Ceo (C6o@CNT) as a func- 
tion of the carbon nanotube (CNT) diameter. For exothermic peapods (CNT 
diameter > 11.8 A), only minor changes, ascribed to a small structural defor- 
mation of the nanotube walls, were observed. These include a small electron 
charge transfer (less than 0.10 electron) from the CNT to the Cgo molecules 
and a poor mixing of the Cqq orbitals with those of the CNT. Decreasing 
the diameter of the nanotube leads to a modest increase of the charge den- 
sity located between the Cgo's. More significant changes are obtained for 
endothermic peapods (CNT diameter < 11.8 A). We observe a large electron 
charge transfer from Cgo to the tube, and a drastic change in electron trans- 
port characteristics and electronic structure. These results are discussed in 
terms of 7r-7r interaction and Cgo symmetry breaking. 
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I. INTRODUCTION 



The possibility of using new forms of carbon-based materials for practical applications 
in nanoelectronics has stimulated an important amount of exciting works during the last 
decade. [1] Since their discovery in 1991, carbon nanotubes (CNTs) have been used as an 
active component in the fabrication of transistor [2-4], memory elements [5], and more re- 
cently, logic circuits. [6,7] CNT's can also be used as a template for nanofabrication and 
as reservoirs for the storage of gas, ions, or metals. [8] In this respect, it has been recently 
shown experimentally that multiple C 60 molecules can penetrate into a carbon nanotube to 
form a one-dimensional array of C 60 nested inside. [9] This new type of carbon materials, 
due to its original structure, is often called carbon peapod. 

Filling CNTs with Cqq is exothermic or endothermic depending on the size of the nan- 
otube. [10] C 60 @(10,10) was found to be stable (exothermic) while other peapods with 
a smaller CNT shell such as the (9,9) and (8,8) tubes are endothermic. Metallic CNTs 
preserve most of their intrinsic properties upon the encapsulation of Cgo molecules. The 
interaction between Cqq and the nanotube occurs through a weak orbital mixing between a 
near free electron (NFE) state on the CNT located above Fermi level (and most probably 
above vacuum level) and the p orbitals of the C 60 . This interaction leads to a weak electron 
confinement between C 60 and the CNT wall, and a slight charge transfer from the tube to 
the Ceo- Therefore, only little perturbation is expected. However, recent STM results show 
drastic modification of the local electronic structure of semiconducting nanotube peapods. 
[11] This perturbation is essentially present in the conduction band. The main variation 
observed is a sharp increase in the density of states probed directly over an encapsulated 
C 60 molecule. This increase in the DOS was tentatively attributed to the electronic cou- 
pling of C 60 with the CNT shell, which is estimated approximately to about 1 eV. [11] This 
tube-Ceo coupling appears much stronger than anticipated by theory [10] or experimental 
work on metallic peapods. [12] It is therefore important to study further the properties of 
encapsulated Ceo on the electrical and electronic properties of both metallic and semicon- 
ducting nanotubes as a function of the nanotube diameter. In the present study, we show 
that the electronic as well as the electron transport properties of exothermic peapods do 
not show drastic differences from the properties of individual species. The CNT-Cgo inter- 
action increases slightly as the nanotube diameter decreases, but becomes very important 
for endothermic peapods. 

II. COMPUTATIONAL DETAILS 

We considered the encapsulation of up to three Ceo molecules in approximately 10 nm 
long metallic and semiconductor CNT models containing up to 1500 carbon atoms. The 
Cgo molecules were systematically placed in the middle sections of the CNT. The region 
where C 60 's are encapsulated is relatively small with respect to the entire CNT length, thus 
avoiding the influence of the open boundary conditions implicit to our finite model. The 
computed stabilization energy (AEs) corresponds to the energy difference between a fully 
optimized peapod structure and the isolated species (perfect CNT and isolated Ceo's). Prior 
to the optimization, the dangling bonds at both ends of the finite nanotube models were 
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saturated with hydrogen. The activation energy needed to introduce the Ceo's in the CNT 
was not evaluated. Stabilization energies were calculated with the MM3 molecular mechanic 
force field [13] (the bond parameter of alkene was modified to 1.42 A) and the geometries 
were optimized with a standard conjugated gradient technique down to a root mean square 
(RMS) deviation < 1CT 5 . We previously showed that this modified-MM3 force field gives a 
total energy that is in good agreement with the more accurate TB-DFT method [14]. The 
electronic structure calculations for the CNT and peapod systems were carried out within 
the extended Hiickel (EH) method, which includes an explicit treatment of overlap integral 
for the s and p valence orbital of carbon. [15] It has been shown that EH gives results similar 
to those obtained on extended CNTs with more sophisticated methods. [16] 

The electrical transport properties of carbon nanotubes and peapods were computed us- 
ing a Green's function approach [17] within the Landauer-Biittiker formalism. The Hamilto- 
nian and overlap matrices used in this formalism were also determined using the EH model. 
[15] For transport calculations, the two ends of the CNT were bonded to gold electrodes. 
The metallic contacts consist of a sufficient number of gold atoms in a (111) crystalline ar- 
rangement to create a larger contact area relative to the CNT ends. In order to minimize the 
contact resistance, the distance between the gold pad and the tube end was fixed to 1.5 A. 
[14] The Green's function of the conductor can be written in the form of block matrices 
separating explicitly the molecular Hamiltonian: 

G c = [ESc-Hc-^i-^T 1 (1) 

where Sc and Tic are the overlap and the Hamiltonian matrices of the conductor (nanotube 
or peapod), respectively, and S 12 are self-energy terms describing the effect of the leads. 
The transmission function T(E) (or transmittance), which is the summation of transmis- 
sion probabilities over all conduction channels in the system, is obtained from the Green's 
function of the conductor (Gc) given by [17]: 

f(E) = T 21 = Tr[r 2 G c I\ G c \ (2) 

In this formula, the matrices have the form: 

r 1>2 = i(E 1>2 - si i2 ) (3) 

III. RESULTS AND DISCUSSION 

The existence of encapsulated C§q inside a CNT can be discussed in term of stabilization 
energy (AE$) in the form of: 

AE S = E[xC 60 @(n,m)\ - E[(n,m)\ - xE[C m ] 

where x is the number of C 60 molecules, and (n, m) is the chirality index of the nanotube 
considered. The variation of the stabilization energy as a function of the nanotube diameter 
(Djvt) for a single encapsulated C 60 molecule is shown in Figure 1, while Table 1 gives AE S 
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(reported by Ceo) values for multiple encapsulated Ceo molecules. In order to qualitatively 
describe the influence of the 7r-electron cloud on the resulting peapods stability, we show an 
additional curve (the dotted line) in Figure 1 where the thickness of the 7r-cloud (~ 3.3 A) is 
substracted from the nanotube diameters. The stabilization energy for peapods with a large 
diameter CNT shell is relatively weak and exothermic. On the other hand, C 60 can more 
easily penetrate into large CNT because the interaction energy is weak. As the CNT diam- 
eter decreases, the interaction energy between the Ceo molecule and the CNT wall becomes 
more important; the AE$ values reflect then the balance between van der Waals attraction 
and Coulomb repulsion. The most stable peapod is the xC 60 @(10, 10) system, in which the 
7r-electron clouds of both C 60 and CNT just begin to overlap. For CNT smaller than the 
(10,10) tube, the peapods become rapidly less stable and highly endothermic (i.e. positive 
AE S ) at DjVT < 11-9 A. The range of tube diameters where peapods are the most stable (12 
< Dnt < 15 A) is in good agreement with recent experimental [18] and theoretical [10,19] 
observations. The C 6 o@(9,9) peapod was previously found [10] slightly endothermic (by 6 
kcal/mol) with DFT-LSD method. The lower stability found for C 60 @(9,9) with DFT is 
probably related to the imposed commensurate structure of the Ceo within the tube in their 
supercell model. The relative stability of small peapods is strongly related to the ability 
of the tube to satisfy the presence of encapsulated C 6 o molecules through a deformation 
of its structure, especially near the C 6 o molecules. As expected, the smaller is the CNT, 
the larger are the structural deformations, and consequently lower is the stability of the 
peapod. [10,19] This also reflects on the Ceo-Ceo distance (dc 60 -c 60 )i C^o m °l ecu les m small 
CNTs become squashed and the distance between their centers increases (see Table 1). In 
the following, we first compare the electronic and electrical properties of the most stable 
metallic (C 60 @(10, 10)) and semiconducting (C 60 @(16,0)) peapods. Then, we present the 
results for the highly endothermic C 60 @(14,0) peapod. 

Figure 2 compares the density of states (DOS) of perfect-CNT/peapods systems (upper 
panels), and the local density of states (LDOS) of the CNT shell of these peapods (lower 
panels) in a region near the C 60 , for the (A) metallic 3C 60 @ (10,10) and (B) semiconducting 
3C 60 @(16, 0) structures. The vertical dotted-dashed lines indicate the energy position where 
a difference in LDOS was observed between a perfect CNT and a peapod for which the 
contribution of the Ceo molecules was removed. This last comparison allows us to highlight 
the changes induced by a structural deformation rather than by an electronic influence of 
C 60 on the CNT. The finite DOS at the Fermi energy (E F ) for C 60 @(10,10), and the absence 
of states at E F for C 60 @(16,0) peapod suggest that the fundamental (metal, semiconductor) 
electronic characteristics of the CNTs are preserved upon C@o encapsulation. The presence 
of Ceo in the peapods clearly results in additional peaks in DOS for valence and conduction 
bands of both metallic and semiconducting peapods. The first two peaks near Ep associated 
to C 60 are observed at around -1.0 eV and +0.2 eV for 3C 60 @(10,10), and at around -0.7 eV 
and +0.5 eV for 3C 60 @(16,0). These peaks correspond respectively to the HOMO and the 
LUMO of C 60 . The calculated HOMO-LUMO gap (1.2 eV) for encapsulated C 60 is smaller 
than the value calculated for isolated Cgo (1-6 eV), and the experimental values (~ 1.6-1.8 
eV) for Ceo in gas and solid phases. [20] This difference is partly related to the deformation 
of Ceo which contributes to closing the gap of Ceo, [22] and to the displacement of Ceo 
orbitals induced by the relatively weak interaction between C 6 o and the CNT. In addition, 
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the very small charge transfer [21] from the CNT to the Ceo molecule (see Table 1), and 
the weak mixing of states between Cqq and the CNT (Cqo states are weakly spread in the 
DOS for peapods) support that only weak interaction between C 60 and the nanotube wall 
are present. This result is in agreement with previous DFT-LSD description of metallic 
peapods, in which a very weak charge transfer was observed. [10] The influence of the CNT 
diameter on the charge density distribution is represented in Figure 3 in which we consider 
the residual charge density p r such as: 

p r = p[xC 60 @(n, to)] - p[xC 60 ] - p[(n, to)] 

and where the dark and bright regions indicate a gain and a loss of electron charges den- 
sity, respectively. As the CNT diameter decreases, the accumulation of charge density 
between the C 60 molecules, and between the C 60 and the CNT shell increases, leaving the 
C 60 molecules slightly more negative than the isolated species. This result contrasts slightly 
with previous DFT result where the weak accumulation of negative charge density in large 
diameter peapods was mainly located in the space between the tube and Ceo- [10] This 
charge density localization, which is practically absent for peapods smaller than the (10,10) 
tube, was attributed to the presence of a weak coupling between C 60 and a near free electron 
(NFE) state of the CNT, which is known to be poorly described within EH method. The 
case of Cgo@(14,0) peapod (D) is quite different. There is an important loss of charge density 
between the Cqo and the CNT shell and a small gain of charge density between Ceo's. As 
discussed below, this behavior is mainly related to Cqo symmetry breaking. 

The LDOS of peapods (lower panels of Figure 2) also suggests a weak influence of the 
Cqo molecules on the electronic properties of the CNT peapods. In these LDOS diagrams, 
the label "0" marks the carbons that are the closest to C6o 5 and the labels ±1, ±2, . . . 
are for carbon sections progressively away from the central "0" position. For the metallic 
3C6o@ (10,10) peapod (see Figure 2A), the most important changes occur in the valence band 
between -1.7 and -1.1 eV. These variations are more directly related to the states created 
by structural deformation of the CNT shell (indicated by the vertical dashed-dotted lines) 
as opposed to the possible electronic influence of Ceo- A similar situation occurs for the 
semiconducting 3Ceo@(16,0) peapod, except that a smaller tube diameter results in slightly 
higher electronic influence on the LDOS of the (16,0) tube. A certain number of structural 
deformed states (dashed-dotted lines) coincides with the presence of C 60 states and suggest 
a possible mixing of states. In contrast to the DOS where the contribution of the tube 
and the three Ceo molecules are convoluted, the local variation of DOS in the vicinity of 
a single Ceo remains almost imperceptible. This result is in strong constrast with previous 
STM measurements in which a very important increase of LDOS was observed but only for 
the conduction band of a semiconducting peapod. [11] The slight variations in valence and 
conduction bands observed in the present study for semiconducting (and metallic) peapods 
suggest that an additionnal phenomenom could occur in the STM experiment and may ex- 
plain the discrepancy. For example, a Coulomb blockade event or a structural deformation 
of carbon peapod induced by the STM tip can alter significantly the spectroscopic signature. 
There is clearly a possible blockade in that the charging energy of C 6 o exceeds 270 meV, 
which is much larger than kT. [23] 
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The transport properties of the (A) metallic (10,10) and (B) semiconducting (16,0) nan- 
otubes are not very much altered by the encapsulation of Cqo's. The main panels of Figure 
4 shows the variation of the transmission function T(E) (or transmittance) for perfect and 
C 60 filled nanotubes as a function of electron energy. We reproduce also the DOS curves of 
CNTs and peapods (lower panels) in Fig. 4 to identify the energy position where changes are 
observed. For the metallic tube, the presence of transmittance peaks, instead of a plateau, 
near Ep is mainly due to our finite model for which the bands still have a molecular (dis- 
crete) character. The imperfect gold-CNT contact induces an extra contact resistance (of ~ 
8 kf2) to the minimal resistance of 6 kf2 for a metallic (ballistic) nanotube, and lowers the 
transmittance at E F from 2 to th 0.8. [24] However, since the CNTs and peapods have similar 
gold-tube geometry, we are expecting a similar contact resistance. The regions where small 
changes in transmittance are observed agree well with the features observed in DOS associ- 
ated to C 60 (see Figure 4A). As observed for the electronic structure properties, the small 
diameter of the (16,0) tube facilitates a slight improvement of orbitals mixing between the 
CNT and the C 60 's. As a result, the transport properties of the (16,0) tube are more altered 
than for the (10,10) tube, especially in the valence band. However, these changes remain 
quite weak and suggest it would be experimentally very difficult to differentiate between 
a perfect nanotube and a peapod, at least on the basis on their electronic and electrical 
properties. 

The existence of the highly endothermic C6o@(14,0) peapod is highly improbable because 
of the large activation energy needed for encapsulation (see Table I). This highly deformed 
peapod is however quite rich in information. The structure of Ceo@(14,0) contains a bumped 
CNT shell near C 60 's. Because of the small inner space in a (14,0) tube, the C 60 molecules 
are significantly squashed into an elliptic shape. Figure 5 shows the variation of the elec- 
tronic (A) and electrical (B) properties between a perfect (14,0) tube and a 3Ceo@(14,0) 
peapod. The DOS and LDOS diagrams (left panels) of the peapod show drastic variations 
with respect to the perfect (14,0) tube, more specifically in the band gap region where 
low energy bound states appear for the peapod. The influence of C 60 can now be clearly 
identified in the LDOS curves over a large range of energy in both valence and conduction 
band. In addition, as the vertical dashed-dotted lines indicate, the variation of electronic 
structure of the (14,0) tube is mostly induced by structural deformations. This is deduced 
from a comparison of the perfect (14,0) tube and the empty peapod system. Considering 
the large charge transfer from C 60 to the tube (see Table 1 and Figure 3D), it is clear that 
C 60 has an important electronic and structural effect on the (14,0) tube properties in a pea- 
pod. This influence reflects on the transport properties where important fluctuations of the 
transmittance are observed near the first band edge of the conduction and valence bands. 
Nevertheless, the change of transmittance at Ep is very weak and only a small difference is 
observed between the perfect tube and the peapod. Although we do not want to emphasize 
the electronic and transport properties of this endothermic peapod, the electronic structure 
of the squashed C 60 in this 3C 60 @(14,0) peapod is useful to gain insight about the nature of 
the charge transfer involved in peapods. 

In Figure 6, we compare the electronic structure of three different geometries of a Ceo 
triad, one with a ideal C 60 geometry (A), a second with the C 60 structure as in the (10,10) 



6 



tube (B), and a third with the compressed Ceo structure as in the (14,0) tube (C). The sym- 
metry group of the orbitals and their corresponding positions in the different arrangements 
are also included. The five-fold degenerate h u and three-fold degenerate t\ u bands corre- 
spond to HOMO and LUMO, respectively. Except for the small displacement of the HOMO 
and LUMO orbitals, and the small symmetry breaking of the g g + h g manifold band at higher 
binding energy, the electronic structure of encapsulated-like Ceo molecules as in (10,10), is 
very similar to free Ceo- On the other hand, the symmetry breaking of Ceo orbitals is very 
significant for the compressed geometry as in the C 6 o@(14,0) peapod. The more important 
band splitting is about the HOMO (h u ) where the symmetry breaking produces one four-fold 
degenerate band, and an isolated orbital that is strongly shifted toward E F . This type of 
band splitting was already predicted for a single Ceo molecule by Joachim and coworkers. 
[22] The presence of an occupied orbital near Ep of the peapod increases significantly the 
ability of C 6 o to donate electrons. In addition, as observed in LDOS of 3C 6 o@(14,0) peapod 
(see Figure 5), the low- lying energy states in the band gap of the (14,0) tube related to 
the regions of deformed CNTs, are placed in the appropriate range of energies to receive an 
extra electron from Ceo- A charge transfer from the low- lying orbital of squashed Ceo near 
Ep to the tube would then explain the important loss of charge density on the Ceo molecules 
observed near the nanotube wall reported in Table 1 and the large net positive charge on 
C 60 for the 3C 60 @(14,0). 

In summary, the electronic and electrical properties of metallic and semiconducting car- 
bon nanotubes for the most stable carbon peapods are not significantly altered by the 
encapsulation of Ceo- Only weak charge transfer are observed from the nanotube wall to 
the C 6 o molecules. This change, mostly located between C 6 o, supports a weak orbital mix- 
ing between the two species. As the nanotube diameter decreases, within the exothermic 
peapods limit, a small increase in charge transfer and orbital mixing is observed. For the 
case of endothermic peapods, the changes in electronic and electrical properties are very 
drastic. The most important and relevant effect remains the Ceo symmetry breaking that 
induces the splitting of the HOMO (h u band) into several components especially one near 
E F . This effect improves the electron donation ability of C 60 in the peapods. 
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Table 1. Variation of the carbon peapods properties as a function of the nanotube 

diameter. 
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FIGURES 



FIG. 1. Stabilization energy AE$ for the introduction of a Ceo molecule into a 100 A long 
carbon nanotube (H-terminated) . Inner and outer diameter is when we considere the thickness 
of the 7r-electron cloud. Exothermic peapods have negative AE$, and endothermic peapods have 
positive values. 

FIG. 2. Density of states (DOS) and local density of states (LDOS) of (A) metallic (10,10) and 
(B) semiconducting (16,0) nanotubes and peapods. DOS (upper panels) compares perfect (dotted 
line) and Cqq filled tubes (full line), while LDOS (lower panels) shows the contribution of a circular 
section of carbon atoms in the vicinity of a Ceo molecule (0 is directly over Ceo, and ±1, ± 2 is 
gradually away from Ceo)- 

FIG. 3. Representation of the residual valence charge density of the 3Cqo@(ti, m) peapods 
where (A) n = m=10, (B) n = 16, m = 0, (C) n = m = 9, and (D) n = 14, m = 0. A negative 
value in the scale corresponds to a loss of charge density (bright region) while a positive indicates 
a gain of charge density (dark region). A similar scale is kept from (A) to (C) to emphasize the 
effect of the tube diameter on the charge density, and the scale for (D) is an order of magnitude 
higher. 

FIG. 4. Variation of the electron transport (upper panels) and electronic (lower panels) prop- 
erties from (n,m) to 3CQo@(n,m) systems, for (A) n = m = 10 and (B) n = 16, m = 0. The 
electronic structure of a perfect (n, m) tube (lower panels) is indicated by the dotted line-streaked 
area. 

FIG. 5. Comparison of electronic structure (A) and electrical (B) properties between the (14, 0) 
and the 3C 6 o@(14,0) systems. 

FIG. 6. Effect of space (tube) diameter on the electronic structure of (A) free Ceo when encap- 
sulated in (B) (10,10), and (B) (14,0) CNT. Vertical lines show the origins and the displacement 
of splitted orbitals induce by symmetry breaking. 
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FIGURE 3. Rochefort 



13 





14 




FIGURE 5. Rochefort 



15 




FIGURE 6. Rochefort 



16 



